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Abstract 

We study the production of W bosons at large transverse momentum in pp collisions at 
the Large Hadron Collider (LHC). We calculate the complete next-to-leading order (NLO) 
corrections to the differential cross section. We find that the NLO corrections provide a 
large increase to the cross section but, surprisingly, do not reduce the scale dependence 
relative to leading order (LO). We also calculate next-to-next-to- leading-order (NNLO) 
soft-gluon corrections and find that, although they are small, they significantly reduce the 
scale dependence thus providing a more stable theoretical prediction. 



1 Introduction 



The study of ly-boson production in hadron colliders is important in testing the Standard 
Model [T] and in estimating backgrounds to Higgs production and new physics 0. Precise 
theoretical predictions for W production at the Large Hadron Collider (LHC) are needed to 
exploit fully the large number of such events that are anticipated when the LHC begins opera- 
tion in 2007. Accurate predictions for this process, which has a particularly clean experimental 
signature when the W decays to leptons, are needed to reduce the systematic uncertainties in 
precision measurements, e.g., of the W mass and decay width. Due to the large number of 
events expected, this process can also be used to determine the parton-parton luminosity jS]. 
Furthermore, accurate calculations for W production at large transverse momentum, Qt, are 
required to distinguish the Standard Model prediction from signals of possible new physics, 
such as new gauge bosons and extra dimensions li4|, which may be expected to enhance the Qt 
distribution at large Qt- 

Calculations of the next-to-leading order (NLO) cross section for W production at large 
transverse momentum at the Fermilab Tevatron collider were presented in Refs. 0ini. These 
predictions have been tested experimentally by the CDF and D0 collaborations [7] , and found 
to be consistent with the data at large Qt- The NLO corrections contribute to enhance the 
differential distributions in Qt of the W boson and they reduce the factorization and renormal- 
ization scale dependence of the cross section. A recent study |HI included soft-gluon corrections 
through next-to-next-to- leading order (NNLO), which provide additional enhancements and a 
further reduction of the scale dependence. 

At leading order (LO) in the Quantum Chromodynamic (QCD) coupling as, a W^-boson can 
be produced at large Qt by recoiling against a single parton which decays into a jet of hadrons. 
The NLO corrections to this cross section involve one-loop parton processes with a virtual gluon, 
and real radiative processes with two partons in the final state. The virtual corrections involve 
ultraviolet divergences which renormalize as and cause it to depend on a renormalization energy 
scale which we will take to be ~Qt- Both real and virtual corrections have soft and collinear 
divergences which arise from the masslessness of the gluons and our approximation of zero 
masses for the quarks. The soft divergences cancel between real and virtual processes. The 
net collinear singularities are factorized in a process-independent manner and absorbed into 
factorization-scale-dependent parton distribution functions. Complete analytic expressions for 
the inclusive NLO cross sections in terms of the partonic Mandelstam variables were presented 
in Refs. and form the basis for the NLO results in this work. 

The calculation of hard-scattering cross sections near partonic threshold, such as W^-boson 
production at large transverse momentum, involves corrections from the emission of soft gluons 
from the partons in the process. At each order in perturbation theory one encounters large 
logarithms that arise from the cancellations between real emission and virtual processes, due 
to the limited phase space available for real gluon emission near partonic threshold. These 
threshold corrections formally exponentiate as a result of the factorization properties jU] of the 
cross section. Expansions of the resummed cross section can be derived in principle to any 
higher order. A unified approach to the calculation of NNLO soft-gluon corrections for hard- 
scattering processes was recently presented in Ref . |10j . Analytic expressions for W production 
are given in Ref. jH|. 
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We note that there has also been work on resummation of Sudakov logarithms at small 
transverse momentum for electroweak boson production in hadron colliders fn\ IT^ IT^ ITl] as 
well as on joint resummation of Qt and threshold logarithms |15j. The results presented in this 
paper for large Qt > 20 GeV are not sensitive to this resummation. 

In this paper we calculate the complete NLO corrections as well as the NNLO soft-gluon 
corrections for VT-boson production at large transverse momentum at the LHC We study the 
size of the corrections and their significance in stabilizing the cross section versus changes in 
the factorization and renormalization scales. 



2 Kinematics and Partonic Cross Sections 

For the production oi a. W boson, with momentum Q, in collisions of two hadrons Ha and hs 
with momenta Pa and Pb 

hA{PA) + hB{PB) W{Q) + X , (2.1) 

where X denotes all additional particles in the final state, we can write the factorized single- 
particle-inclusive cross section as 



X Eq ^^I^^^^is, t, u, Q^ /iF, asif^D) (2.2) 

where Eq = Q^, the parton distribution (pf/h is the probability density in the momentum 
fraction x for finding parton / with momentum p = xP in hadron h, and a is the perturbative 
parton-level cross section. The initial-state coUinear singularities are factorized into the parton 
distributions at factorization scale ^p, while /i/j is the renormalization scale. 

At the parton level, the lowest-order subprocesses for the production oi a,W boson and a 
parton involve quarks anti- quarks q and gluons g: 

qiPa) + g{Pb) WiQ) + qipc), 

q{Pa) + q{Pb) W{Q)+g{p,). (2.3) 

The partonic kinematical invariants in the process are s = {pa+PbY, t = {pa—QY,u = (pb—Q)'^, 
and S2 = s + t + u — Q'^ = {pa + Pb — QY ■ Here S2 is the invariant mass of the system recoiling 
against the electroweak boson at the parton level and parameterizes the inelasticity of the 
parton scattering, taking the value S2 = for one-parton production. 
Schematically, the parton level cross sections at NLO have the form 

Eq '^^ = S{s2)as{ijji) [A{s, t, u) + as{iJR)B{s, t, u, /i/j)] + a^(/i^)C(s, t, u, Ss, /if) • 

(2.4) 

The coefficient functions A, B and C depend on the parton flavors fafb and on electroweak 
parameters (which are suppressed), in addition to the kinematic variables s,t,u, S2, and scales 
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^R, fiF shown explicitly. The coefficient A{s,t,u) arises from the LO processes in Eq. ()2.3|) . 
The functions B and C result from a perturbative QCD calculation of the virtual and real 
corrections, respectively, to the processes ()2.3|) . The NLO calculation is done using dimen- 
sional regularization and modified minimal subtraction (MS) to deal with ultraviolet, soft, and 
coUinear divergences. The function B{s,t,u, fipt) is the finite (after renormalization) sum of 
virtual corrections and of singular terms ~(5(s2) in the real radiative corrections. The function 
C{s, t, u, S2, ^f) is the finite (after factorization) contribution from real emission processes away 
from the edge S2 = of phase space. Analytic expressions for these functions from Ref. |6j were 
used to compute the numerical NLO results in this paper. 

In general, the partonic cross section a includes distributions with respect to S2 at n-th 
order in the QCD coupling of the type 



ln™(s2/g|) 



S2 



m < 2n — 1 



(2.5) 



defined by their integral with any smooth function / by 
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dS2f{s2) 
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Qt 




[2.6) 



These "plus" distributions are the remnants of cancellations between real and virtual contribu- 
tions to the cross section. Below we will make use of the terminology that at n-th order in 
the leading logarithms (LL) are those with m = 2n — 1 in Eq. ()2.5|) . the next-to- leading loga- 
rithms (NLL) are those with m = 2n — 2, the next-to-next-to-leading logarithms (NNLL) are 
those with m = 2n — 3, and the next-to-next-to-next-to-leading logarithms (NNNLL) are those 
with m = 2n — 4. Below the term "NNLO-NNNLL" means that the soft-gluon contributions 
through NNNLL to the NNLO corrections have been included and added to the complete NLO 
result (for details see Ref. [S]). 



3 Numerical results 

We now apply our results to W production at large transverse momentum at the LHC Through- 
out we use the MRST2002 parton densities jTH]. The QCD coupling as(/w|./A^) is computed 
using the appropriate value of the 4-fiavor QCD scale parameter A4 required by the LO, NLO, 
and NNLO parton densities, and assuming five fiavors of quarks (u,d,s,c,b). The coupling 
is evolved using the massless LO, NLO and NNLO /^-function as appropriate, and using the 
standard P effective-fiavor scheme to cross heavy-quark thresholds. The electroweak coupling 
a(M|) is evaluated at the mass of the Z boson, and standard values yyj are used for the various 
electroweak parameters. 

In all of the numerical results we present the sum of cross sections for W~ and produc- 
tion. W bosons at the LHC will be detected primarily through their leptonic decay products 
e.g., W~ — > iug. Our cross sections should therefore be multiplied by the appropriate branching 
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pp-->W S =14TeV |Li=Q^ 




Q^(GeV) 

Figure 1: The differential cross section, da ldQ\^ for W production in collisions at the LHC 
with 75 = 14 TeV and = ^ip = = Q^. Shown are the LO, NLO, and NNLO-NNNLL 
results. 

ratios (~0.11 for each lepton species). Our cross sections are integrated over the full phase space 
of the W and final state partons. Detailed comparisons with experiment will require realistic 
cuts on lepton and jet momenta. These modifications will affect the overall normalization of 
the cross sections, but they should not affect the main results of this paper concerning the 
convergence and reliability of perturbative corrections. 

In Fig. [T]we plot the transverse momentum distribution, da/dQ'^, for W production at the 
LHC with "s/S* = 14 TeV. Here we set /ij? = /ir = Qt and denote this common scale by /x. We 
plot LO, NLO, and NNLO-NNNLL results. In the LO result we use LO parton densities, in 
the NLO result we use NLO parton densities, and in the NNLO-NNNLL result we use NNLO 
parton densities. We see that the NLO corrections provide a significant enhancement of the LO 
Qt distribution. The NNLO-NNNLL corrections provide a further rather small enhancement 
of the Qt distribution, which is hardly visible in the plot. Part of the reason for the small 
difference between the NLO and the NNLO-NNNLL curves is that the NNLO parton densities 
reduce the cross section. As we will see below, the NNLO-NNNLL corrections can be much 
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pp-->W S =14TeV |i=Q, 
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Figure 2: The i^-factors for the differential cross section, da/dQ\, for W production in pp 
collisions at the LHC with = 14 TeV and fi = fip = jJ'R = Qt- Shown are the K- factors 
for NLO/LO and NNLO-NNNLL/NLO. The latter is shown with the same or different parton 
distribution functions (pdf) for the NLO and NNLO results. 

bigger for other choices of factorization and renormalization scales. We also note that the 
NNLO-NNNLL result is practically the same as at NNLO-NNLL (i.e. if we had kept one level 
less of logarithms), and that while the LL, NLL, and NNLL terms are complete, the NNNLL 
terms are not complete but include the dominant contributions (see |S] for details). 

In Fig. 121 we plot the i^'-factors, i.e. the ratios of cross sections at various orders and 
accuracies to the LO cross section, all with fip = fiR = Qt, in the high-Qy region. We see that 
the NLO/LO ratio is rather large; the NLO corrections increase the LO result by about 30% to 
50% in the Qt range shown. In contrast, the NNLO-NNNLL/NLO ratio for this choice of scale 
is rather small. Part of the reason for this is that the NNLO parton densities are significantly 
smaller than the NLO parton densities. To make this point more clearly we also show a curve for 
the NNLO-NNNLL/NLO K-factor where both the NNLO-NNNLL and the NLO cross sections 
are calculated using the same NNLO parton densities. The ratio is then significantly bigger. 
We can also see that the NNLO-NNNLL/NLO iiT-factors are nearly constant over the Qt range 
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Figure 3: The differential cross section, da/dQ^, for W production in pp collisions at the 
LHC with y/S = U TeV, Qt = 50 GeV, and /i = /i^ = fiR. Shown are the LO, NLO, and 
NNLO-NNNLL resuhs. 



shown even though the differential cross sections themselves span three orders of magnitude in 
this range. 

In Figs. El in andElwe plot the scale dependence of the differential cross section for Qt = 50, 
80, and 150 GeV, respectively. We set fip = and denote this common scale by /i. We plot 
the differential cross section versus h/Qt over two orders of magnitude. We note that the scale 
dependence of the cross section is not reduced when the NLO corrections are included, but 
we have an improvement when the NNLO-NNNLL corrections are added. The NNLO-NNNLL 
result approaches the scale independence expected of a truly physical cross section. 

It is interesting, and perhaps a little surprising, that the scale dependence of the LO results 
in Figs. EllllElis less pronounced than that of the NLO results. One might expect that the scale 
dependence should decrease with each successive order of perturbation theory. This is indeed 
what is observed at Tevatron energies [HIIH], where the LO cross section increases monotonically 
with decreasing scale and has positive curvature. The LO curves in Figs. 01 El El have little 
curvature, and at Qt = 80 GeV the curve actually turns over at small h/Qt- In Fig. El we 
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Figure 4: The differential cross section, da/dQ^, for W production in pp collisions at the 
LHC with y/S = U TeV, Qt = 80 GeV, and /i = /i^ = fiR. Shown are the LO, NLO, and 
NNLO-NNNLL resuhs. 



plot the LO scale dependence separately for jjp and fi^ with the other held fixed. The cross 
section increases with positive curvature as the renormalization scale is decreased: this is 
the expected behavior at LO due to asymptotic freedom. The dependence however has 
negative curvature and the cross section increases with scale. This behavior is due to the fact 
that the cross section is dominated by the gluon-initiated process qg Wq. The gluon density 
in the proton increases rapidly with scale at fixed x smaller than ~0.01. At LHC energies, the 
fiR and fip dependencies cancel one another approximately, while the scale dependence of the 
LO cross section at Tevatron energies is dominated by the fiR dependence of as- 

In Fig. [7|we plot the differential cross section da/dQ"^ at high Qt with \^ = 14 TeV for 
two values of scale, Qt/2 and 2Qt, often used to display the uncertainty due to scale variation. 
We note that while the variation of the LO cross section is significant and the variation at NLO 
is similar to LO, at NNLO-NNNLL it is very small. In fact the two NNLO-NNNLL curves lie 
very close to or on top of each other. These results are consistent with Figs. IH El 
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Figure 5: The differential cross section, da ldQ\^ for W production in collisions at the LHC 
with \/5 = 14 TeV, Qt = 150 GeV, and fi = fip = l^R- Shown are the LO, NLO, and 
NNLO-NNNLL resuhs. 



4 Conclusions 

We have presented the full NLO and the NNLO soft-gluon corrections for W production at large 
transverse momentum in pp collisions at the LHC. We have shown that the NLO corrections are 
large but do not diminish the scale dependence of the cross section relative to LO. The NNLO- 
NNNLL corrections are very small for /x = Qt but they significantly decrease the factorization 
and renormalization scale dependence of the transverse momentum distributions. These precise 
and reliable perturbative QCD predictions for inclusive W production will facilitate precision 
tests of the Standard Model, measurements of parton densities, and searches for new physics 
at the LHC. 
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Figure 6: The differential cross section, da/dQ^, for W production in pp collisions at the LHC 
with y/S = U TeV, Qt = 80 GeV. Shown are the LO results with fip = /^_r, and with fip and 
fiR varied separately. 
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